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Propagation of Fast Partons in the Nuclear Medium 
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Abstract. The color dipole approach has been applied in the target rest frame to address the issues of 
transverse momentum broadening and energy loss of a fast quark propagating in the nuclear medium. A 
recent application of the theory to the FermiLab E772/E866 experimental data, determining the rate of 
energy loss of a quark propagating in the medium to be 2 to 3 GeV/fm, will be reviewed. Calculations 
for the transverse momentum distribution will be presented, and the results will be compared to the E866 
data. The theory will be shown to compare favorably to the data, and these results will be shown to suggest 
that the momentum broadening of a quark is about twice the generally accepted size. 

PACS. 13.85.Qk; 24.85+p Drell-Yan processes; medium effects; heavy-ion collisions 



1 Introduction 



Nucleon-nucleus (or deuteron-nucleus) reactions at high 
energy are important experimentally because they consti- 
tute the conventional background for the less well known 
dynamics of relativistic nucleus- nucleus collisions, which 
have been identified as a means for producing the quark- 
gluon plasma in the labratory at facilities such as the Rela- 
tivestic Heavy Ion Collider (RHIC) at BNL, and the Large 
. <^ _ Hadron CoUider (LHC) at CERN. To develop sufficient 
, confidence in theory to interpret the more complicated 
5^ ' phenomena in nucleus-nucleus collisions, it is vital that 
pi I the theory be tested against these simpler data first. 

Shadowing, quark energy loss, and transverse momen- 
tum broadening are three central issues of parton propa- 
gation in nucleus-nucleus collisions that may be quantita- 
tively addressed in p+A reactions, specifically in Drell-Yan 
(DY) reactions. The color dipole approach in the target 
rest frame provides an attractive means for interpreting 
these data to learn more about these phenomena. After a 
brief description of the status of this theory, I will review 
recent results that are improving our understanding of the 
physics behind these aspects of parton propagation in the 
nuclear medium. 

Drell-Yan in the target rest frame may be described as 
a quark of the incident hadron interacting with a target 
nucleon, radiating a virtual photon 7* of mass M that 
subsequently decays into the observed Drell-Yan dilepton 
pair In the target rest frame, the 7* is a constituent 
of projectile Auctions, which are "frozen" by time dilation 
for a length of time tc given by the uncertainty relation. 



tc 



2E„ 



(1) 



where Eq and niq refer to the energy and mass of the 
projectile quark and M'^j^ is the square of the effective 
mass of the fluctuation, 
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where a is the fraction of the light-cone momentum of the 
incident quark carried by the lepton pair, and k"^ is the 
square of the transverse momentum of the lepton pair. The 
fluctuation lifetime is called the coherence time, and 
the DY reaction occurs when the 7* of the fluctuation is 
released by an interaction between one of the constituents 
of the fluctuation and a target nucleon. In the color dipole 
approach, this interaction is mediated by the color dipole 
cross section, Oqq. 

When the scattering takes place on a nucleus, multi- 
ple interactions with target nucleons can give rise to var- 
ious medium effects. Two limiting cases should be distin- 
guished, the short coherence length limit (SCL) reached 
when the coherence length ^c = c{tc) is much smaller than 
the interparticle spacing £c <S^ d (in a heavy nucleus, 
d ~ 2fm), and the long coherence length limit (LCL) 
reached when ic ^ Ra, where Ra is the nuclear radius. 
The coherence length determines the physics and also con- 
trols the way that medium effects are taken into account. 
The theory simplifies these important limits. The interme- 
diate case is generally more difficult to describe; however, 
if < ^ Ra, the result can be obtained by interpolating 
between these limits using the square of the longitudinal 
form factor, F^{qc), where Qc = 1/lc is the longitudinal 
momentum transferred in the reaction. The Green func- 
tion method ^ was developed to handle the more difficult 
situations. Our results make use of all these cases. 

One important medium effect is shadowing, and a ba- 
sic consideration is the connection between shadowing and 
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the coherence length. In the SCL there is no shadowing be- 
cause the duration of the fluctuation is so short that it has 
no time to interact with the medium. In the LCL, which 
is applicable to reactions at the LHC and at RHIC under 
certain kinematic conditions, there is maximal shadowing. 

Medium effects also give rise to quark energy loss and 
to transverse momentum broadening. These both arise in 
the SCL as the quark from the incident hadron under- 
goes various additional interactions with nucleons of the 
nucleus before it radiates the 7*. For the LCL, the quark 
may undergo interactions that lead not only to shadowing 
but also to additional momentum broadening. Nuclear ef- 
fects giving rise to momentum broadening and shadowing 
are believed to arise predominantly from the same color 
dipole cross section <Jqq{rT) that mediates the DY reaction 
on a nucleon. 

The DY data relevant to shadowing, energy loss, and 
momentum broadening in which we are interested here 
is the p + A FermiLab data from the E772/E866 collab- 
oration. This data corresponds to the range of coherent 
length where the shadowing is weak and the interpolation 
formula described above applies. 



2 Drell-Yan Reaction on a Nucleon in the 
Target Rest Frame 

In the color dipole approach, the Drell-Yan reaction oc- 
curs when a fast projectile quark scatters off the gluonic 
field of the target, as shown in Fig. ^ The color dipole 
approach was originally proposed for deep-inelastic scat- 
tering (DIS). Extension to Drell-Yan was developed by 
Kopeliovich 2 and subsequently by Brodsky, et al. (Sj- 
The DY reaction on a nucleon is described by 

(3) 

where Uqq is the color-dipole cross section, p is the trans- 
verse distance between the 7* and quark in the fluctu- 
ation, and represents the distribution for the incident 
quark to fluctuate into a quark and the 7*. In DY, the 
color dipole consists of the quark before and after the re- 
lease of the 7*, whose impact parameters differ by ap. 
Here Xq{xi) is the fraction of light-cone momentum of 
the incoming hadron h (for us a proton) carried by the 
quark (lepton pair); the light-cone momentum fraction a 
is given in terms of these variables, a = xi/xq] and, F^^ is 
the quark distribution function of the incident hadron. 

Deep-inelastic scattering from nucleons at HERA at 
DESY has been used to fix models of the color-dipole 
cross section. The HERA data suggests saturation, the 
condition that the color dipole cross section approaches a 
constant value at large r^, 

a,-,(rT) = ao(l-e-'-'/«°). (4) 

Models that incorporate this saturation property include 
the GW W and KST 'S^ models. We see that in both the 




Fig. 1. In the target rest frame, DY dilepton produc- 
tion looks like bremsstrahlung. A quark or an anti-quark 
from the projectile hadron scatters off the target color 
field (denoted by the shaded circles) and radiates a mas- 
sive photon, which subsequently decays into the lepton 
pair. The photon decay is not shown. The photon can be 
radiated before or after the quark scatters. 



GW and KST models, for small tt, 

(^qqirr) « Cr^, (5) 

where C = ctq /i?Q is a constant that depends on the quark 
energy, in the case of the KST model, and on Bjorken x 
in the case of the GW model. For us, the most important 
point is that once the color dipole cross section has been 
fit to experimental DIS data on the nucleon, the color- 
dipole formalism makes predictions for Drell-Yan both in 
nucleon-nucleon W and nucleon-nucleus scattering with 
no further adjustment of parameters. 

The theory is expected to be valid for small Bjorken x 
only, and specific calculations [H] have shown it to agree 
well for X2 < 0.1 in both the magnitude and shape with 
the next-to-leading order parton model. Although the cal- 
culations disagreed with the E772 FermiLab data, they are 
in excellent agreement with the more recent E866 analysis 
inCl. 

3 Medium Effects in DY Reactions for p+A 
Collisions in the Target Rest Frame 

Medium effects are evident in nuclear ratios of cross sec- 
tions R-^/^ . These include analyses of both cross sections 
themselves, and the analysis of momentum distributions. 
For each, there are two cases of interest: weak shadowing 
and maximal shadowing. 

As we have stated, shadowing is controled by the co- 
herence length, and explicit calculations in [5] show that 
for incident quarks of the 800 GeV projectile protons the 
coherence length is quite a bit less than the nuclear radius, 
implying weak shadowing. Additionally, these calculations 
show that the coherence length of the dominant q—^* fluc- 
tuation is nearly constant as a function of xi and is there- 
fore essentially a function of X2 alone. For X2 > .03, ic < 2 
fm, and shadowing can be safely ignored. For smaller X2 
the coherence length exceeds the interparticle spacing at 
the center of the nucleus, and shadowing therefore begins 
to become important. The particular admixture of long- 
and short-coherence length contributions is determined by 
the longitudinal form factor, as stated in the introduction, 
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and this depends in turn explicitly on the variation of 
with xi and X2. 

For cross section ratios, 
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DY 



dAPdxi ' dAPdxi 



(6) 



the E772/E866 FermiLab data at Ep = 800 GeV has been 
used for a recent determination of quark energy loss [5]. 
The DY reaction at RHIC for a;i?=a;i — a:2>0.5 and 
at the LHC lies in the regime of long-coherence length 
physics, and such experiments at these facilities will be 
exploring a limit that has so far not been accessible ex- 
perimentally. Predictions of cross sections inp + A, d + A, 
and A + A collisions have been made recently in the 
latter case. 

For momentum distribution ratios, 



'DY 



'DY 



dp^ 



dp^ 



(7) 



data relevant to the the case of weak shadowing has been 
taken at Ep = 800 GeV by the FermiLab E772/E866 col- 
laboration and was used by them to study quark momen- 
tum broadening in nuclei, 



(8) 



This same regime is relevant io p + A collisions at RHIC. 
The case of maximal shadowing is relevant to the LHC 
and RHIC, and predictions are made in Ref. 0. 



3.1 Cross Section Ratios 

Let us consider first the case of weak shadowing for the 
cross section ratios, as analyzed in |H]. To calculate the 
cross section, we must calculate the cross section sep- 
arately in the SCL and LCL and then interpolate be- 
tween them using using the longitudinal form factor, as 
explained in the introduction. 

The main medium effect in the short coherence length 
limit X2 3> 0.03 the 7* is the quark energy loss, which 
occurs before the 7* is radiated. The interactions that 
cause the incident quark to lose energy reduce the quark 
momentum fraction Xq and hence shift the value of xi = 
aXq. The amount of reduction due to quark energy loss 



may be calculated as Xq « Eq/Eh 



KL/Eh, where 



L is the path length of the incident quark in the nucleus 
and K is its rate of energy loss. Using these kinematic 
considerations, we find an expression in |H] for the DY 

cross section dAF^'dxi hicluding its dependance on k. 
It becomes sensitive to k for xi ~ 1. 

Mechanisms of energy loss jS] consist of the primary 
ones of string breaking (SB), for which « 1 GeV/ fm 
(numerically, the string tension) , and the independent pro- 
cess of gluon radiation (GR), for which kgr ~ 3as{k'^) « 
0.8 GeV/ fm, where w 0.65 GeV is the mean-square 
momentum of the radiated gluon. The induced energy-loss 



mechanisms arising from multiple interactions with nucle- 
ons are relatively minor for nuclei of ordinary density and 
can be neglected. We thus expect that k « 1.8 GeV/ fm. 

As we have said, shadowing arises in the LCL through 
multiple interactions between the quark and the nucleus 
following the emission of the 7*. For weak shadowing, 
Glauber theory to second order in the number of inter- 
actions may be used, giving 



dM^dxi 



^eff{TA) ■ 



(9) 



where the brackets around the color dipole cross section 
refer to the averages over (p, a) in Eq. and where (Ta) 
is the average of the thickness function TA{b) over impact 
parameter b, 



< Ta > -- 



A 



dHTlib), 



with TA{b) = TA{b, 00), where 

Ta(6, z) = / d'^bpA{b,z), 



(10) 



(11) 



and pA{b,z) is the nuclear density. The quantity (Jeff — 
{(Jqq) / {(Jqq) in Eq. jnj is obtained by expanding the mul- 
tiple scattering through second order in aqq. In the color- 
dipole approach shadowing is completely determined by 
the color dipole cross section when shadowing is weak. 

Interpolating between the LCL where X2 ^ 0.03 and 
the SCL where X2 ^ 0.03, we obtain 



d ody o" 



.2 (SCL) 



DY 



dIvPdxi dRPdxi 



■(1-Fl(g,)) + 



2 (LCL) 



DY 



dAPdxi 



(12) 

where Qc = 1/ic and the longitudinal form factor is defined 
as 



Fliq) 



1 



<Ta> 



d'^b\ / dze''''' p{b, z)\' 



(13) 



Note that this interpolation formula has the correct limits 
for £c ^ and £c ^ 00. It also correctly reproduces the 
Gribov shadowing formula 10.1.11) and correctly describes 
shadowing in DIS ■S-- Because the cross section becomes 
sensitive to both shadowing and k for a;i « 1, it is impor- 
tant that both effects are included in the theory (note that 
the depencence of shadowing on £c becomes quite intricate 
in this limit 12,8 ). This theory then provides a quan- 
titative test of the value of the relatively poorly known 
quark energy loss in the nuclear medium, since shadowing 
is completely determined by theory. 

For comparing the theory to experiment, we have eval- 
uated Eq.((nil and compared to the E866/E772 FermiLab 
data. This data provides the dependence of R^/^ over a 
range of {A,xi, M^); our ability to reproduce these data 
gives us confidence that we are able to separate the effects 
of shadowing and energy loss. Typical results are shown 
in Fig. El 
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Fig. 2. Examples of ratios of the DY cross sections of 
tungsten to deuterium as functions of xi for various inter- 
vals of M. Dashed curves correspond to net shadowing 
contribution, solid curves show the full effect including 
shadowing and energy loss. Data are from Refs. \lHll4f . 



The rate of energy loss corresponding to the fit to the 
the data, as shown in Fig. |21 is 



2.73 ± 0.37GeV/fm 



(14) 



Note that the value in Eq. H14|) is somewhat larger than 
the theoretical estimate of 1.8 GeV/fm given above; it 
is also considerably larger than values found in previous 
analyses of the DY data ■ We regard our value as more 
reliable, since the color dipole approach is able to disentan- 
gle shadowing from energy loss using a reliable theoretical 
calculation of the shadowing contribution. 



3.2 Momentum Distribution Ratios |16 | 

Calculations show that the momentum distributions from 
E772/E866 experiments with SOOGeV' protons or measure- 
ments at RHIC with s^/^ = 200 GeV protons at a;F = 
correspond to the short coherence length limit, as long as 
X2 > 0.05. Working in the short coherence length limit 
simplifies the theory by eliminating the shadowing term 
in Eq. 1)12(1 . We will restrict our attention to this kine- 
matic regime, which is possible because the data is avail- 
able binned in X2 intervals, and thus concentrate on the 
first (LCL) term in this expression. 

It was recognized in Ref. J7] that at high energies, the 
multiple interaction of a quark propagating in the nuclear 
medium can be eikonalized and exponentiates into a factor 
containing the color dipole cross section and nuclear thick- 
ness function. The physics of momentum broadening of the 
quark is therefore that of color filtering, or absorption of 



large dipoles leading to diminished transverse separation 
with distance. Then, the probability distribution W^{kT) 
that a quark will have acquired transverse momentum fcr 
at a position (6, z) in the nucleus A becomes [T7) 



1 



d ttg 



ikT-TT „-^<yqq{rT)TA(b,z) 



(15) 



where Ta(&, z) is defined in Eq. Hll|) . We express the cross 
section, a^yiPT): for a- proton to produce a DY pair on a 
nucleus A with transverse momentum px as convolution 
of the probability W^{kT) and the DY momentum dis- 
tribution on a proton. The resulting expression entails an 
integral over !Z'(a, p) as in Eq. Q, and for the results we 
show below, this integral has been performed numerically. 
However, for the purpose of discussing the result here, we 
replace the quantities that depend on a by averages. Then, 
the ratio R^/p{pt) may be expressed as in terms of the 
DY cross section crDYiPr), 



'DY 



(pt) 



J d^kTW'^ikT)alY{PT - akr) 



where a arises from the fact that the DY pair carries away 
a fraction a of the transverse momentum of a quark in 
the incident proton. We find that a — (a^)^/^ « 0.97. 
The average over a implicit in Eq. H16|) also entails an 
evaluation of the quark energy, (Eq) « Eq/3, since we use 
the GW model of aqq, which depends explicitly on this 
energy. 

For our calculations below we have chosen 



aP,YiPT)^ii+PT/y^y 



(17) 



where the value of = 7GeV^ is taken from Ref. |18| . 

Some of the integrals in Eas. (|15ll6|l may be done an- 
alytically by expanding out the exponential in Eq. 
We find that to an excellent approximation, the effect of 
averaging Eq. (fTK|l over (6, z) is to replace 



TA{b,z) ^ {Ta)/2. 



(18) 



When evaluating Eq. (|15|l . we also take into account gluon 
shadowing by replacing aqq RaCTqq, where Rq is the 
gluon shadowing function. See Ref. for more discus- 
sion. Our calculations are however insensitive to Rq- 

Our calculation for R^^^'^{pt) is shown in Fig. El The 
data relevant to pt distributions on nuclei are FermiLab 
E772/E866 data [TMniil9, . One sees that the calculation 
is in quite good agreement with experiment. Note that the 
details of the Cronin effect, the rise above 1 for ~ 3 in 
Fig. 13 are completely explained with no adjustable param- 
eters. We find comparable agreement with R^^/^'^{pt)- 



3.3 Transverse Momentum Broadening 

The transverse momentum broadening of a quark propa- 
gating in a nucleus is defined in Eq. |(SJ) where the mean 
momentum is defined as 



{pD 



f dlj,p^j.a^Y{PT) 



pT'^DY 



(pt) 



(19) 
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Fig. 3. Comparison of theoretical prediction of 
R^''^'^{Pt) vs. pt (in GeV/c) to experiment for X2 ~ 
0.05. Data are from the FermiLab E772/E866 collabora- 
tion li,'^IJ4li.9t . 



The dependence of the quantity 6{p'^) on A has been stud- 
ied experimentally using the E772 data in Ref . . They 
found that for large A, 



T/exp 



0.021A^/^GeV^. 



(20) 



We can make a calculation of the same quantity in the 
theory just described. From Eq. (|16|l we easily find that 



<^dy{pt) = — 



d\Te'P--^W{rT)aly{rT) 



where 
and 



e 2 



\cri,{rT)TA{b,z) 



'DY 



(21) 



(22) 



(23) 



characterizes how a quark acquires transverse momentum 
as it propagates through the nucleus. From Eq. p9l21|l we 
find 



2 \A 



72 „P 



{Pt) 



'DY 



>t)| 



<^DYi^T) 



V^W^^(arT) 



\rT=0 (24) 



The first term in Eq. |(23Jl, the mean-square momentum 
on a nucleon, is formally infinite. Fortunately, the nucleon 
contribution gets subtracted to obtain Sip'^), 



Hpt) 



\rT=0- 



(25) 



Note that the shadowing contribution in Eq. (|12|l . which 
we justified in dropping by going to X2 > 0.05, leads to a 



divergent contribution that no longer cancels exactly; in 
cases where FSI are important other means are needed to 
define a meaningful quantity 9 . 

Using Eq. H22|l . along with the observation that crqq{rT) 
is proportional to r|, at small (see Eq. (|5}), 



5{pI) = a^C{TA) 



(26) 



where we have used the average of as given in Eq. H18|) . 

We may now make estimates to compare with the ex- 
perimental result in Eq. (|20|l . Using values of C w 4 — 5 
as obtained from the GW or KST models, and taking the 
sharp-surface model for the nuclear density. 



{Ta) = -^PoRa 



(27) 



where po w 0.16/m ^ is the central density of heavy nuclei 
and Ra ~ l.lAa, we find from theory that 



S{p't) 



Th 



0.047^3 GeV^ 



(28) 



The theoretical result in Eq (|28|l is clearly larger than and 
in apparent disagreement with the experimental result in 
Eq. 1201 17,21 . However, it has been pointed out |221 
that large systematic errors must be included in Eq. (|2()|l 
and that when these are taken into account, the result in 
Eq. could be accommodated. 



3.4 Predictions for Maximal Shadowing 

As we have stated, the case of maximal shadowing cor- 
responds to ic 3> Ra, and experiments in this regime 
have not yet been possible. The color dipole approach in 
the target rest frame applies in this limit, and numerous 
predictions have been presented in Ref. 1^. We give one 
example in Fig. ^ namely the DY ratios as a function 
of transverse momentum presented separately for longitu- 
dinal (L) and transverse -I- longitudinal {T + L) photon 
polarization. 

In contrast to the weak shadowing limit discussed ear- 
lier, transverse momentum broadening here is associated 
with interactions with the quark following the emission of 
7*. Note that the Cronin peak is much less visible in the 
LCL, which is a consequence of the gluon shadowing. For 
evaluating the gluon shadowing here, the Green's function 
method is required, and we take Rg from the calculations 
of Ref. EBl. 



4 Summary and Conclusions 

We have shown that the color dipole approach formulated 
in the target rest frame has advantages for examining nu- 
clear modifications in p+A collisions because the effects of 
shadowing may be determined theoretically. In the SCL, 
£c Ra, energy loss and momentum broadening are 
mediated by multiple interactions of the incident quark 
with target nucleons before the emission of the 7*. Analy- 
sis of the E772/E866 FermiLab experimental DY data at 
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Fig. 4. Nuclear effects on the DY transverse momentum 
distribution. Curves show the DY cross sections for pAu 
(left) and deuterium - gold (right) collisions per nucleon 
divided by the DY cross section from pp scattering. Solid 
curves are predictions for RHIC (y'^ = 200 GeV) and 
dashed for LHC (^ = 5.5 TeV). Calculations are for 
M=4.5 GeV andxF= 0.5. 



Ep = 800 GeV in the color dipole approach has provided 
improved determinations of these quantities, giving 



dEq 
dz 



-2.73 ±0.35 GeV/fm 



and 



(5(p|} « 0.047yl3 GeV"^ 



(29) 



(30) 



The value of the rate of energy loss is larger than previous 
determinations, and the value of S{p'^) is about twice as 
large as the conventional value. 

Predictions of the transverse momentum dependence 
of DY cross sections, for p + A and A + A collisions, have 
been made in the color dipole approach. These provide 
tests of the theory in a new physical regime of strong shad- 
owing, where £c ^ Ra, that may be explored at the LHC 
in the future. 
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